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Abstract. A number of classical algorithms are based on random walks on 
graphs. It is hoped that recently defined quantum walks can serve as the basis 
for quantum algorithms that will faster than the corresponding classical ones. 
We discuss a particular kind of quantum walk on a general graph. We affix two 
semi-infinite lines to a general finite graph, which we call tails. On the tails, the 
particle making the walk simply advances one unit at each time step, so that its 
behavior there is analogous to free propatation. We are interested in how many 
steps it will take the particle, starting on one tail and propagating through the 
graph (where its propagation is not free), to emerge onto the other tail. The 
probability to make such a walk in n steps and the hitting time for such a walk 
can be expressed in terms of the transmission amplitude for the graph, which 
is one element of its S matrix. Demonstrating this neccessitates a study of the 
analyticity properties of the transmission and reflection coefficients of a graph. 
We show that a graph can have bound states that cannot be accessed by a 
particle entering the graph from one of the tails. Time-reversal invariance of 
a quantum walk is defined and used to show that the transmission amplitudes 
for the particle entering the graph from different directions are the same if the 
walk is time-reversal invariant. 



1. Introduction 

Random walks on graphs serve as the basis of algorithms for solving a number 
of problems, including 2-SAT, graph connectivity, and finding satsifying assign- 
ments for Boolean functions. Quantum algorithms have shown promise in solving 
some problems faster than is possible using classical algorithms. Quantum walks 
represent an attempt to "quantize" classical random walk algorithms and thereby 
increase their speed. 

Quantum algorithms must be run on a quantum computer. In this kind of 
machine, information is represented not in terms of bits, but in terms of quantum 
bits, or qubits. A bit is either or 1 and can be represented by current on or 
off, or by two different voltages. A qubit is a two-level quantum system, and 
possibilites include the spin states of an electron or the polarization states of a 
photon. One of the states, e.g. spin up, corresponds to and the other, e.g. spin 
down, corresponds to 1. The principle difference between a bit and a qubit, is that 
the bit is either or 1, while the qubit can be in a superposition of its state and 
its 1 state. Computers based on bits represent and process information according 
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to the rules of classical physics, while computers based on qubits represent and 
process information according to the rules of quantum physics. 

There has been some progress in using quantum walks to construct quantum 
algorithms. It has recently been shown that it is possible to use a quantum walk to 
perform a search on the hypercube faster than can be done classically pQ. In this 
problem the number of steps drops from N, which is the number of vertices, in the 
classical case to y/~N in the quantum case. A much more dramatic improvement has 
recently been obtained by Childs, et al. 2 . They constructed an oracle problem 
that can be solved by a quantum algorithm based on a quantum walk exponentially 
faster than is possible with any classical algorithm. 

Quantum walks were introduced by Y. Aharanov, L. Davidovich and N. Zagury 

They considered a spin-1/2 particle that is shifted to the left or right depending 
on its spin state. Since then, a number of different kinds of quantum walks have 
been developed. The first, the coined walk, is a discrete time walk that makes use 
of an auxiliary quantum system, the coin. This is necessary to make the operator 
that advances the walk one step unitary. This walk was originally proposed by 
Watrous and is closlely related to the original walk discussed by Aharanov 0|. A 
second type of walk is the continuous time walk, due to Childs, et al. 0. Here the 
Hamiltonian that is used to construct the time development transformation is taken 
to be proportional to the adjancency matrix of the graph. A third type of walk, 
based on thinking about the graph as an interferometer with optical multiports as 
the nodes was recenty proposed by us 0. In this case the walk takes place on 
the edges of the graph, rather than the vertices, and each edge has two states, one 
corresponding to traversing the edge in one direction and the other to traversing 
the edge in the opposite direction. It should be noted that while it is very simple 
to construct a quantum walk on an arbitrary graph using the second and third 
methods, it is much less obvious how to do so for the coined walk. 

Let us briefly describe the coined quantum walk. In trying to formulate a 
quantum walk on a graph, the most natural thing to do is to let a set of orthonormal 
basis states correspond to the vertices of the graph. If a particle is in the state |n), 
that corresponds to its being located on vertex n. Let the Hilbert space spanned by 
these states be Tip. Trying to define a unitary evolution using this scheme soon leads 
to serious problems, as was first noted by Meyer [Zj. Watrous solved this problem 
by enlarging the Hilbert space in which the quantum walk takes place. How this 
scheme works is most easily seen by considering the quantum walk on a line. The 
vertices are labelled by integers, and, in addition, there is a quantum coin, which 
has two states, \L) and \R), corresponding to left and right, respectively, which 
span a two-dimensional Hilbert space, Ti c . A basis for the Hilbert space describing 
the entire system, Ti p (8 H. c , is given by the states \n) ® \a), where n is an integer, 
and a is either L or R. A step in this walk consists of applying the Hadamard 
operator, H, to the coin, 
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where S is the shift operator, whose action is given by 
(1.3) S|n) = |n + 1) &\n) = \n-l). 

Note that the evolution of the quantum state of the walk is deterministic, that is, 
if the initial state is |^, n ), then the state after m steps is just [Vh(I p <8> H)] m \ip in ). 
Probabilities enter the picture when we try to determine where the particle is by 
making measurements. 

Coined quantum walks have been studied on both the line by Nayak and Vish- 
wananth |8] and on the cycle by D. Aharanov et al. 9 . Aharonov et al. also 
studied a number of properties of quantum walks on general graphs. Absorbing 
times and probabilities for walks on the line have been examined by several authors 
1 1 01 111) . Going beyond one dimension, properties of quantum walks on two and 
three dimensional lattices |12) and on the hypercube I14| have been explored. 
The effects of decoherence on quantum walks have also been studied. Brun, et al. 
showed how increasing amounts of decoherence turn a quantum walk into a classical 
random walk |15j . Kendon and Tregenna found that small amounts of decoherence 
can speed up the convergence of the time-averaged probability distribution of a 
quantum walk |16| . Many aspects of quantum walks are covered in the recent and 
excellent review by Kempe |17| . 

Several proposals have been made for the physical realization of quantum walks 
1181 |24| . The ones closest in spirit to the walks studied here are the realizations 
that employ optical methods, either linear optical elements |21L I22j or cavities 
[231 124) . The last two references show that an experimental quantum walk has, in 
fact, been carried out, though it was not intrepreted as such at the time |25| . 

Here we wish to describe quantm walks that take place on the edges of a graph, 
and to study the connection between these walks and scattering theory. In order 
to show how scattering theory enters the picture, begin with a graph, choose two 
vertices, and attach two semi-infinite lines, which we shall call tails, to it, one to 
each of the selected vertices. Each of the half-lines is made up of vertices and edges. 
Thus we can start the walk on one of the tails, have it progress into the original 
graph, and emerge onto the opposite tail. This type of arrangement allows us to 
define a scattering matrix, or S matrix, for the original graph, with the amplitude 
to get from one tail to the other being called the transmission amplitude. As we 
shall see properties of quantum walks starting on one tail and ending on the other 
can be expressed in terms of the transmission amplitude of the graph. It should be 
noted that scattering approaches have also proven useful in a related area, quantum 
graphs, which can be used to study electron transport in large molecules |26| . 

The first use of scattering theory to analyze the behavior of quantum walks was 
done by Farhi and Gutmann |27| . They studied continuous-time walks on trees, 
and they added tails to the trees to place the walks into a scattering-theoretic 
framework. They were able to place bounds on the time for a walk, starting at the 
root of the tree, to penetrate to one of the leaves. 

We begin by defining our quantum walk, and then illustrate it with an example. 
We then explore the relation between properties of the quantum walk and the 
transmission amplitude of the graph. These are then applied to an example. We 
define a time-reversal transformation for quantum walks and use it to explore when 
the transmission amplitude to go through the graph in one direction is the same 
as the transmission amplitude to go through it in the other direction. Finally, we 
finish with a detailed study of the properties of the transmission and reflection 
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amplitudes, in particular showing that they are analytic in on the unit circle and 
in its interior. This allows us to express the probability of starting on one tail and 
finishing on the other in n steps and also the hitting time for a walk from one tail 
to the other very simply in terms of the transmission amplitude of the graph. 

2. Quantum walks on edges 

The type of quantum walk that will be used throughout this paper was origi- 
nally presented in [5] . We imagine a particle on an edge of a graph; it is this particle 
that will make the walk. Each edge has two states, one going in one direction, the 
other going in the other direction, and we denote the set of oriented edges (an edge 
and its direction of traversal is an oriented edge) of the graph by E. That is, if our 
edge is between the vertices A and B, which we shall denote as (A, B), it has two 
orthogonal states, \A, B), corresponding to the particle being on (A, B) and going 
from A to B, and \B,A), corresponding to the particle being on (A, B) and going 
from B to A. The collection of all of these edge states is a basis for a Hilbert space, 
and the states of the particle making the walk lie in this space, which is L 2 (E). 

Now that we have our state space, we need a unitary operator that advances 
the walk one step. Let us first illustrate how this works for a walk on the line. 
We shall label the vertices by the integers. In this case, the states of the system 
are \j,k), where k = j ± 1. The vertices can be thought of as scattering centers. 
Consider what happens when a particle, moving in one dimension, hits a scattering 
center. It has a certain amplitude to continue in the direction it was going, i.e. to be 
transmitted, and an amplitude to be change its direction, i.e. to be reflected. The 
scatterer has two input states, the particle can enter from either the right or the 
left, and two output states, the particle can leave heading either right or left. The 
scattering center defines a unitary transformation between the input and output 
states. 

We now need to translate this into transition rules for our quantum walk. 
Suppose we are in the state \j — 1, j). If the particle is transmitted it will be in the 
state \ j, j + 1) , and if reflected in the state \ j, j — 1) . Let the transmission amplitude 
be t, and the reflection amplitude be r. We then have the transition rule 

(2-1) \j-hj)^t\j,j + l)+r\j,j-l), 

where unitarity implies that \t\ 2 + \r\ 2 = 1. The other possibility is that the particle 
is incident on vertex j from the right, that is it is in the state \j + If it is 

transmitted it is in the state — and if it is reflected, it is in the state | j, j + 1) . 
Unitarity of the scattering transformation then gives us that 

(2.2) \j + l,j)^t*\j,j-l)-r*\j,j + l). 

These rules specify our walk. 

The case t = 1 and r = corresponds to free particle propagation; a particle 
in the state \j, j + 1) simply moves one step to the right with each time step in the 
walk. If r ^ 0, then there is some amplitude to move both to the right and to the 
left. As an example, let us consider the case when t = r = 1/V2, and the particle 
starts in the state |0, 1). The probability distribution for finding the particle on a 
given edge after 1, 000 steps is shown in Figure 1. First, note that this distribution 
look nothing like the distribution that would result from a classical random walk 
on a line, which would be a Gaussian centered about the initial position of the 
particle. In addition, we see that there is an interval about the origin in which it is 
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Figure 1. Probability distribution for quantum walk after 1000 steps 

very likely to find the particle. An asymptotic analysis of the walk shows that this 
region lies between — \t\n and \t\n, where n is the number of steps in the walk and 
t is the transmission amplitude of the vertices ^] . That is, the region in which the 
the particle is likely to be found grows as n instead of \fn as would be the case in 
a standard random walk. 

Thinking about the particle making the walk and the scattering centers as beam 
splitters, we see that this type of quantum walk is analogous to the propagation of 
light through an interferometer. This suggests that we can add a new element to 
quantum walks that has no analogue in classical random walks. Interferometers are 
made up of multiports (generalized beam splitters) j2g] and phase shifters; a phase 
shifter imparts a constant phase to a photon that passes through it. Suppose we 
were to put a phase shifter that imparts a phase shift of (j> just before the j th vertex. 
The transition rules for the states adjacent to this vertex are modified, while the 
rules for all other states are unaffected. In particular, we now have 



Insertion of a phase shifter into an edge can change the properties of a quantum 
walk, because it changes how different paths interfere. 

So far we have only considered vertices at which two edges meet, but if we are 
to construct graphs more complicated than lines, we need to see how a vertex with 
more that two edges emanating from it behaves. We shall begin by looking at an 
example. 

If a vertex treats all edges entering it in an equivalent fashion, then we have 
a particulary simple situation, because the edges of the graph do not have to be 
labelled. One vertex of this type is very closely related to the quantum coin used in 
the walk on the hypercube |13l I14 | . Let the vertex at which all of the edges meet 
be labelled by O, and the opposite ends of the edges be labelled by the numbers 
1 through n. For any input state, |/c,0), where k is an integer between 1 and n, 
the transition rule is that the amplitude to go the output state \0,k) is r, and 
the amplitude to go to any other output state is t. That is, the amplitude to be 
reflected is r, and the amplitude to be transmitted through any of the other edges 



(2.3) 



b'-i,j) 
b' + i,j) 



te^\j,j + l)+re 2 ^\j,j-l) 
-r*\j,j + l)+t*e^\j,3-l). 
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is t. Unitarity places two conditions on these amplitudes 

(n-l)|i| 2 + M 2 = l 

(2.4) (n - 2)\t\ 2 + r*t + t*r = 0. 

We shall call such vertices equal-transmission vertices. As an example, for the case 
n = 3, possible values of r and t are r = — 1/3 and t = 2/3. 

In order to construct a walk for a general graph, one chooses a unitary operator 
for each vertex, i.e. one that maps the states coming into a vertex to states leaving 
the same vertex. More precisely, let ta be the set of oriented edges starting at the 
vertex A and Ta be the linear span of ta ■ Similarly, let to a be the set of oriented 
edges ending at the vertex A, and Ha be the linear span of lua- We have that 
T A Pi t b = an d UA{~\u B = $iiA^B, which implies that 

(2.5) E = \Jt a = \Ju;a, 

A A 

and 

(2.6) L*(E) = @T A = ®Sl A , 

A A 

where the unions and direct sums are over all vertices of the graph. The local 
unitary operator, Ua maps Ha to Ta, and describes the scattering at vertex A. 
One step of the walk consists of the combined effect of all of these operations; the 
overall unitary operator, U, that advances the walk one step is constructed from the 
local operators for each vertex. Explicitly, the edge state \A, B), which is the state 
for the particle going from vertex A to vertex B, will go to the state Ub\A, B) after 
one step, where Ub is the operator corresponding to vertex B. This prescription 
guarantees that the overall operation is unitary, in particular, U acting on any 
other edge state \C,D) will give a state orthogonal to U\A,B). If D = B, then 
\A, B) and \C,D) will be mapped into Tb, but the unitarity of Ub ensures that 
U\A, B) and U\C, D) are orthogonal. IfD^B, then U maps \A, B) and |C, D) onto 
different sets of states, Tb and To, respectively, and the results are then orthogonal. 
Therefore, as the edge states make up an orthonormal basis of the Hilbcrt space in 
which the walk occurs, and U maps this basis to another orthonormal basis, it is 
unitary. 

3. Example 

Let us put all of this together in a very simple example. Consider the graph 
shown in Figure 2, where each of the vertices where two edges meet have t = 1 and 
r = 0, while the three-edge vertices are of the type discussed in the previous section, 
with r = —1/3 and t = 2/3. There is a phase shifter with an adjustable phase, 
<p, in one of the edges. The graph goes to negative infinity on the left and plus 
infinity on the right. Note that it is very simple to construct the unitary operator 
that advances a quantum walk on this graph by one step; one simply combines the 
actions of the operators corresponding to each vertex. 

We can find the unnormalized eigenstates for this graph, and one set of them 
can be described as having an incoming wave from the left, an outgoing transmitted 
wave going to the right, and a reflected wave going to the left. A second set will 
have an incoming wave from the right, an outgoing trasmitted wave to the left, and 
a refelected wave to the right. The mathematical status of these objects will be 
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Figure 2 . Graph for example 



explored in Sections 7 and 8. Finally, there may be bound states, i.e. eigenstates 
that are localized in the region between the two vertices with three edges. Now, let 
us denote the left three-edge vertex by and the right one by 2, and number the 
vertices on the lines correspondingly, from 2 to plus infinity to the right and from 
to minus infinity to the left. The eigenstates with a wave incident from the left 
take the form 

-l 

l*> - E (e^|i,i + l) + &-ie- i(j+1)e |i + l,j)) + |*o 2 ) 

j=-oo 

CO 

(3.1) +E«2e^'" 2 ^|i,i + 1), 

where |\&02) is the part of the eigenfunction between vertices and 2, and e~ lB 
is the eigenvalue of the operator U that advances the walk one step. The first 
term can be thought of as the incoming wave; it is confined to the region between 
negative infinity and 0, and consists of states in which the particle is moving to the 
right. The term proportional to b-i is the reflected wave. It is also confined to the 
region between negative infinity and 0, but consists of states in which the particle 
is moving to the left. Finally, the term proportional to a 2 is the transmitted, being 
confined to the region from 2 to infinity, and consisting of states with the particle 
moving to the right. The coefficient b-i can be interpretted as the amplitude for 
the incoming wave to be reflected, and a-i as the amplitude to be transmitted. 
Denoting the upper vertex between and 2 as \A and the lower one as IB, we can 
express ^02} as 

|* 02 ) = a 0A \0,lA)+a 1A \lA,2) + b 0A \lA,0) + b 1A \2,lA) 

(3.2) +a OB \0, IB) + a 1B \lB, 2) + b 0B \lB, 0) + b 1B \2, IB). 

The solution of the equation Z7|\E') = e is given in the appendix. 

One can define a transmission coefficient for this graph, which is just the ratio 
of the intensities of the outgoing transmitted and the incoming waves, 

T _l„,2_ 4(l + e-'*)(l-e-W*>) 

{ ' 121 e - 4 ^(l + e-^) 2 -(3e^( 4e +«-l) 2 ' 

One finds that for <f> = 0, T is nonzero, because the two paths from to 2 interfere 
constructively, while when (f> — it they interfere destructively, which results in 
T — 0. Therefore, the behavior of the walk strongly depends on the value of the 
phase shift. 
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4. Hitting times 

Suppose we want to find out whether the particle making the walk is on the 
edge between vertices j and j + 1. The relevant projection operator is 

(4.1) Pj = \j,j + j + 1| + \j + + 1, j|. 

If we obtain 1 the particle is on that edge, and if we obtain 0, it is not. What we 
wish to find out is, if the particle is initially in the state ^(O)), the probability 
that the particle is not on the edge for the first n steps of the walk, but is on the 
(n + l)st. We measure Pj after every step, and we have to take the effect of these 
meaurements into account in describing the dynamics of the wave function. 

Let us first see what happens after one step. Let U be the unitary operator 
that advances the walk one step. The probability that the particle is not on the 
edge between j and j + 1 after one step, p{l), is 

(4.2) p(l) = \\(I-P j )U*(0)f. 

and the quantum state, assuming the particle was not found on this edge, is 

(4.3) |*(i)) = iLi3)Emi. 

- Pj)U*(0)\\ 

Now let us go one more step to see the pattern. The probability that the particle 
is not on the edge between j and j + 1 after either the first or second steps, p(2), is 

P (2) = p(i)\\(i-Pj)u^(i)f 

(4.4) = \\(I-P 3 )U(I-P 3 )U^(0)\\ 2 , 

and the quantum state, assuming the particle was not found on the edge after cither 
step, is 

(i-Pj)umi)) 



(4.5) |tt(2)) = 



\(I-Pj)U*(l) 



In general, if p(n) is the probability of not finding the particle on the edge between 
j and j + 1 after each of n steps, then 

(4.6) p{n) =p(n- 1)||(7 -Pj)U^{n- 1)|| 2 , 

where — 1)) is the quantum state after n — 1 steps and n — 1 measurements 
indicating that the particle was not on the designated edge. We also have that 



(4.7) |*(n)) 



(I-PjpMn-l)) 



\(I -Pj)U^(n- 1)|| ' 
Let us now proceed by induction. We shall assume that 

p(n) = }\[(I-P j )U] n nO)\\ 2 
[(I-P 3 )U} n \^(0)} 



(4.8) |tf(n)) = 



i[(/-p f )tn n *(o)ir 



which clearly holds for n — 1,2. Substitution of the second of these equations into 
Eq. l|4~7Jl yields 



(4.9) |*(n + l)) = 



{I-P 3 )U\^(n)) 
||(/-P,0^*(n)|| 
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so that part of our induction hypothesis is verfied. From Eq. I|4.6J) we find that 

P (n+i) = m-PAurmfw-PAvvwf 

(4.10) = m-p^ur+^m 2 , 

which proves the remaining part. This result is easily generalized to the case of 
making measurements to determine whether the particle is on a set of edges instead 
of on a single edge. A similar result for the coined quantum walk was derived by 
Bach, et at 

The hitting time for a random walk on a graph is the expected number of steps 
in a walk that starts at one specified vertex and stops on first reaching a second 
specified vertex. In order to calculate this we need the probability that a walk 
starting on the first vertex does not reach the second for its first n — 1 steps, but 
does reach it on the nth step. The corresponding quantity for a quantum walk, 
which we shall call q(n), is the probability that a walk starting in the state |^(0)) 
is not measured to be on the edge between j and j + 1 for the first n — 1 steps, but 
is measured to be on that edge at the nth step. This probability is given by 

q(n) = p(n- l)||P.,-Z7#(n - l)f 

(4.11) = \\P j U[(I-P j )U] n - 1 y(0)\\ 2 . 

The particle making this walk may never be found to be on the edge between j and 
j + 1 at all, and the probability that it is, is given by 

oo 

(4.12) Pout = Y, ( ln- 

71=1 

We can define a conditional hitting time for this walk, h, as 

1 °° 

(4.13) h=—J2 nc ln- 

n— 1 

This quantity provides an an indication of how many steps a walk that does reach 
the edge between j and j + 1 will require. 

Let us now consider the situation depicted in Figure 3. We start the walk on 
the edge | — 1, 0) moving to the right, and we are interested in how long it takes to 
get to the edge |j, j + 1), where G is a general graph with a finite number of edges 
and vertices. The vertices outside of the graph, G, for example, —1 and j + 1, have 
a transmission amplitude, t = 1. We shall denote the entire graph, G plus the tails, 
as G. The case considered in the previous section is a special case of this situation. 
When calculating the probability, q(n), that the walk arrives at the edge |j, j + 1) 
for the first time after n steps, we need to evaluate (see Eq. (14.1111 ) 

n-l 

(4.14) PjU[(I - PjW]"- 1 = P 3 U n - PjU( n - k) PjU k + ... 

We first note that once the particle leaves the graph, G, it moves one unit to the 
right with each step. That means that terms of the form PjU m PjU k \ — 1,0), for 
m > 0, are zero, because the first Pj puts the particle in the state \ j, j + 1), and the 
operator U m moves it m steps to the right, so that the particle in now in the state 
\j + m, j + m + 1). The projection Pj acting on this state gives zero. Therefore, 
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Figure 3. A general graph, G, with tails attached to two if its vertices 

only the first term on the right-hand side of Eq. (|4.14() gives a nonzero contribution. 
That means that in this case, we have 

(4.15) q(n) = (*(0)|((7t)«P J [/«|M/(0)). 

We now want to relate the probability q(n) to the transmission coefficient of 
the graph, which we shall define shortly. The eigenstates of the system can be 
grouped into three sets, those with the particle coming in from the left, those with 
the particle coming in from the right, and bound states. The eigenstates of this 
system corresponding to the particle coming in from the left are 

-l 

= (e l{k+1)0 \k,k + l)+r(e)e^ k+1 ^\k + l,k)) + \^ G (9)) 

k— — oo 
oo 

(4.16) +Y / t{0)e l{k -i )e \k,k + l), 

k=j 

where r(0) is the reflection coefficient, t(9) is the transmission coefficient, and 
\^g(9)) is the part of the eigenfunction inside the graph G. This state is, in fact, 
not in L 2 (E), and its mathematical significance will be clarified in a later section. 
In Appendix B we show that 

(4.17) |r(0)| a + |t(0)| 2 = l. 

For the initial state of our walk we will choose | — 1, 0), and we want to express 
this state in terms of the eigenstates in Eq. H4.16fl . To do so we integrate \^i(9)) 
from to 27r, finding 

r 2TT 1 r 2-ir 



±-l d£\M0)) = 1-1,0) + ^/ dff[ £ r(0) e -«* +1 > s |fc + l,fc>] 

(4.18) +— / de[\^ G (9))+Y,t(0y {k - j)e \k,k+l)}. 

71 Jo k=j 

By setting z = exp(«6*) the integrals on the right-hand side of the above equation 
can be turned into integrals along the unit circle, which we shall denote by C, in 
the complex plane. These integrals are of the form 

dzz m - Y r{z) J c dzm G {z)) 

c 



(4.19) / dzz m ^t(z) 

c 
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where m > 0. As we shall show in Section 7, r(z), t(z), and l^ciz)) are analytic 
functions of z on the unit circle and in its interior, and they all vanish when z = 0. 
This implies that all of the above integrals are zero. Therefore, we can conclude 
that 

f-271 



(4.20) 

Now we have that 
(4.21) 



1,0) 



1 
2~/7 



d0|*z(0)>. 



q(n) = \(j,j + l\U n \-l,0)\ 2 , 



and 



(j,j + l\U n \-l,0) 



1 

2^ 



dde- in9 (j,j + l\M0)) 



(4.22) 

Therefore, we find 
(4.23) 



2 71 



j- I d0e~ mS t(9). 



q(n) 



2-it 

— I dde- tr " 
2tt 



t{9) 



which relates a walk probability, the probability to land on the edge between j and 
j + 1 for the first time after n steps, to the S matrix of the graph G. 

We can rewrite the integral appearing in the above expression for q(n) as a 
contour integral over the unit circle, C by setting z — e 10 . This gives us 



(4.24) 



1 

2^ 



d6e 



-inO 



t{9) 



1 

2ni 



1 



c 



dz rrUz). 



Because, as we shall see, t{z) is analytic in a region containing the unit circle and 
its interior, the only contribution to the integral in Eq. 14.24JI comes from the pole 
at the origin that is due to the l/z n+1 factor. This fact makes it straightforward 
to evaluate the above integral, by just finding the residue of the pole at the origin, 
giving us 



(4.25) 



1 

27ri 



c 



dz —rt(z) 



1 d n 
n\ dz n 



t(z) 



z=0 



This is just the n th coefficient in the Taylor expansion of t(z) about the origin, i.e. 
if 



(4.26) 



t(z) 



n=0 



then we have that q(n) = |c„| 2 . The conditional hitting time, defined in Eq. I|4.13ll . 
is then given by 

oo 

(4.27) h = Y,n\c n \ 2 P ov 

±OUt r, 



We now note that 
(4.28) 



*/ s dt 

t{z)z Tz 



oo oo 



rn — n— 1 



nc m Cn(z ) z , 
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from which it follows, setting z = cxp(z#) and integrating, that 
Similarly, we have that 

(4.30) Pout = ^- I** d6\t(e ie )\ 2 . 

The integrals in the previous two equations can also be expressed as contour inte- 
grals if we define the "reflected" transmission amplitude, tn(z) to be tn{z) = t(z*)*~ . 
We then have that 

p °^=hl c dz (S) tR ($) t{z) 

( 4 ' 31 ) h 2iriP out j c dztR (z) dz 

Thus, we see that both the probability of first hitting the edge between j and j + 1 
in n steps and the hitting time to reach this edge are related to the transmission 
amplitude of the graph. 

5. Application to example 

Let us now apply this formula to a walk on the graph discussed in Section 3. 
The transmission amplitude for this graph is 

4e- J0 (l + e-^)(l - e -*(40+«) 



e-™"(l + e-^) 2 - (3e- 4 ( 49 +« - l) 2 ' 



(5-1) M e - m 

and 

(b ' t[Z) ~ z 4 (l + e-^f - (3e-** - z 4 ) 2 " 

From this expression we see explicitly that t{z) is analytic inside and on the unit 
circle. 

Considering the case (j) = we find the simple expression 

Hz 3 

(5-3) t(z) = — 4 . 

Inserting this result into the expressions in the previous section, it is straightforward 
to perform the integrations, and we find that 

g („) = ( (gT^Tiw) 2 ifn = 3mod4 
\ otherwise 

(5.4) P out = | 

For small n, it is relatively straightforward to verify the expression for q(n) by 
adding up the amplitudes for the possible paths from | — 1, 0) to |2, 3) of length n. 

Another issue that deserves mention is the existence of bound states. In stan- 
dard quantum scattering theory, if a potential has bound states, these are orthgonal 
to all scattering states. What we mean by a bound state of a graph, is an eigenstate 
of U whose support is confined to G. If our graph has bound states, then these 
states could not be accessed by a particle coming in from one of the tails, e.g. from 
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the left. This follows from the fact that the initial state is orthogonal to any bound 
state (a state localized on one of the tails is automatically orthogonal to one that is 
localized in G), and if ^i n u is the initial state, and ^bnd is the localized eigenstate 
with eigenvalue exp(— iOf,), we have 

(5.5) (*bnd\U n * init ) = e- ine »(* bnd \* init ) = 0. 

This limits the walks that particles with initial states of this type can perform, 
because, as the above equation shows, the state of the walk at any time must be 
orthogonal to all of the bound states. 

In order to find the bound states in our example we consider a state \^f) that 
is located on the four edges connecting the vertices 0, 1A, IB, and 2, i.e. 

|*) = aoA|0,L4) +a 1A \lA,2) +b OA \lA,0) +b 1A \2,lA) 

(5.6) +aos|0, IB) + a 1B \lB, 2) + b OB \lB,0) + b 1B \2, IB). 

We then solve the equation U\^) — e |\&). This equation only has a solution 
when = 0, and in that case the eigenvalues are i m , where m = 0, . . . 3. The 
eigenstates in this case are given by 

1 

o-oa 



aiA = 

(5.7) a 1B = 

A simple basis for the space spanned by these eigenstates is 



2V2 


b A = 


_ »"• 
2V2 


1 


b\A = ■ 


_ ( _ 1)m 


2i m V2 


2V2 


-1 






2V2 


boB = 


i m 
2V2 


-1 


biB = 


( _ 1)m 


2i m V~2 


2V2 





l«i> = 


^(|0,1A>- 


|0,1B)) 




\U2) = 


i(HA2) - 


|1S,2» 




\ u 3) = 


^(|1B,0>- 


11-4,0)) 


(5.8) 


\ Ui ) = 




|2,l-4» 



A state resulting from a walk starting on the tails, will be orthogonal to the subspace 
spanned by these four states. A walk starting in this subspace will never leave the 
subspace, and hence, will never escape into the tails. 

For all of the eigenvalues except those given by i m , the corresponding space 
of (unnormalized) eigenvectors is two dimensional. It is spanned by two scattering 
states, one with a incoming wave from the left, and one with an incoming wave 
from the right. These states are orthogonal. At the eigenvalues i m the space of 
eigenvectors becomes three dimensional. In addition to the two scattering states, 
there is one bound state. All of these state are mutually orthogonal. 
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6. Time reversal 

In quantum one-dimensional scattering problems, the time reversal invariance 
of the Hamiltonian is used to prove that the transmission amplitude for a wave 
coming from the left is the same as that for a wave coming from the right. Here we 
would like to see if something similar is possible for quantum walks. 

Our first task is to define a time-reversal operator. The standard time-reversal 
operator is anti-unitary, does not affect positions, and reverses velocities. This 
suggests that we define the anti-unitary operator, T, where the hat is used to 
indicate that the operator is anti-linear, by 

(6.1) f \A,B) = \B,A). 

Note that T is its own inverse. For a walk to be time-reversal invariant, we want 
the operator that advances the walk one step, U , to be taken into its inverse when 
conjugated by T, 

(6.2) fUT = U- 1 = U\ 

that is, conjugation with T changes the operator that moves the walk forward one 
step into the one that moves it backwards one step. 

Let us see what this condition implies about U. Let k enumerate the vertices 
that are connected to A by an edge (one of which is B). This set will be denoted 
by T(A). We find that 

TUf\A,B) = ^2(A,k\U\B,A)"\k,A) 

k 

(6.3) U'^B) = J2{A,B\U\k,A)*\k,A) 

k 

These will be the same if 

(6.4) (A, k\U\B, A) = (A, B\U\k, A), 

so that U will satisfy the condition of time-reversal invariance if for all vertices A, 
and k,k' € T(A), we have 

(6.5) (A, k\U\k', A) = (A, k'\U\k, A). 

We find that equal-transmission vertices satisfy this condition, and they still satisfy 
it if phase shifters are put into the edges. We shall henceforth confine our attention 
to graphs that do satisfy this condition. 

If |\E') is an eigenstate of U with eigenvalue exp(— i9), and U is time-reversal 
invariant, then 

(6.6) UT\ty = e- l6 f\y), 

where we have made use of Eq. (|6.2|) . This implies that if |\E') is an eigenstate of 
U, then so is T\^), with the same eigenvalue. Note that if \^f) is a bound state, 
this equation implies either that is invariant under the action of T, or that the 
eigenvalue exp(— iQ) is degenerate. 

We can now show that the transmission coefficients from both directions are 
equal. Consider two eigenstates, each with eigenvalue exp(— iO), one incident from 
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the left, and one incident from the right 

-l 

= E (e 4(fc+1)e |A : ,fc+l>+r i (e) e -^ +1 " ? |fc + l,A : )) + |* Gi (0)) 



k— — oo 
oo 



k—j 

-1 

l*r(0)> = J2 U(0)e^ k+1)9 \k + l,k) + \^ Gr ) 



k— — oo 

OO 

(6.7) ^( e - l(fe - j)e |fc + 1, fc) + r r (9)e^ k - j ^\k, k + 1)). 
fc=j 

Now, consider the eigenstate 

(6.8) |*) = A|*,)+B|* r ). 

If we apply T to this state, we must obtain an eigenstate of U with an eigenvalue 
of exp(— iff), which must, therefore, be a linear combination of |\&/) and I^V), 

(6.9) f|*)=C|*i)+D|* r ). 
This implies that 

A* = Cri + Dt r A*tf + B*r* = D 

(6.10) A*rf + B*t* r = C B* = Cti + Dr r . 

Substituting from the second and third of these equations into the first, gives us 

(6.11) \rt\ 2 +t r tf = l nt*+t r r* =0. 

The first of these equations, and the fact that |r/| 2 + \ti\ 2 = 1, gives us that t r = tj. 
Therefore, time-reversal invariance does imply that the left and right transmission 
amplitudes are the same. 

7. Analyticity of the reflection and transmission amplitudes 

We now want to see how to find, schematically, the reflection and transmission 
amplitudes of a general graph, and then deduce some of their properties. In par- 
ticular we want to show that they are analytic in a region that contains the unit 
circle and its interior. 

We begin by considering the graph G and denote the set of oriented edges of G 
by E. The space in which the walk takes place is Tt = L 2 (E), and this immediately 
presents us with a problem, because \^i(9)) (see Eq. (|4.16|I % 1. through which the 
reflection and transmission amplitudes are defined, is not in L 2 (E). It can certainly 
be viewed as a complex- valued function on the set E, (A, B) € E — > (A, B\^i(6)), 
and it is in this sense that we take the following calculations until we have more 
information on the nature of the specific functions. 

The transmission and reflection amplitudes are found by solving the eigenvalue 
equation, 

(7.1) U\%) = exp(-i0)|tf). 
A direct calculation shows that this is equivalent to 

(7.2) U(\ -1,0) + |* G (0))) = e- w (\^ G (e)) + r(0)|O, -1) + t(fi)\j,3 + !))• 
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Now let Hg be the subspace of H spanned by the oriented edges of G, H3 be the 
subspace of H. spanned by the oriented edges {| — 1,0), |0, —1), 1)}, Pg be the 
orthogonal projection onto the space TIg, and P3 be the orthogonal projection onto 
the space H3. In addition, define H = 7^3© 7Yg and P as the orthogonal projection 
onto Ti. We can view the elements of TL as column vectors of complex numbers 

/ vi \ 

V3 
V v G J 

where v± is the coefficient of | — 1, 0), V2 is the coefficient of |0, —1), V3 the coefficient 
°f I Ji J + 1) an d v<3 the 2iV-tuple of coefficients of the edges of G (we assume that 
G has TV edges). In this notation, Eq. (|7.2|) becomes 



/ 



(7.3) 



where 



(7.4) 



-i0 





r(9) 
t(9) 

VG 



\ 



X 

w 



Y 
G 



V VG J 



( 1 \ 




V V G J 



X 

w 



Y 
G 



= {PUP) = 



P3UP3 
P G UP 3 



PsUPg 
PgUPg 



In the above equation, X is a 3 x 3 matrix, Y is a 2N x 3 matrix, W is a 3 x 2N 
matrix, and G is a 2N x 2N matrix. The second and third columns of both X 
and W are identically zero. We also note that the first rows of both X and Y are 
identically zero. 

It is possible to break the solution of Eq. (|7.3(l into two parts. Denote the single 
nonzero column of W (the first one) by w, i.e. 

/ 



(7.5) 



w = W 



V J 

where is a 27V-tuple of zeroes, so that w G Hg- 
is solving the equation in TCg, 



The first step in solving Eq. H7.3(l 



(7.6) 



(e" ie I - G)v G = w, 



where I is the 2N x 2N identity matrix. Once we have found from this equation, 
we insert the solution into Eq. (|7.3() in order to find r(8) and t(9). 

We can extend the reflection and transmission amplitudes into the complex 
plane by replacing exp(— i6) in the above equations by 1/z. This will yield the 
In particular we have that if vg(z) is a solution to 

(I — zG)vg = zw, 



functions i(z) and r(z) 
(7.7) 



then 



(7.8) 



t(z) 
r(z) 



zY U,j+iy,(k,l) V G;(k,l)(z) 

(k,l)eG 



— Z [^(0 : -1);(-1,0) + 



(kA)eG 



Y(0-l);(k,l)VG;(k,l)( z )}i 
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where the sums are over all edges in G. 

We are first interested in determining whether t(z) is analytic inside the unit 
circle. This will be the case if Eq. (|7.6|) can be solved for all \z\ < 1, which will be 
true if the matrix I — zG is invertible for \z\ < 1. This follows from the result that if 
an operator, T(z), is analytic in some region of the complex plane, so is its inverse, 
T~ 1 (z) provided it exists 29 . Because G = P G UP G we have that ||G|| < 1, which 
implies that if \z\ < 1, then I — zG is invertible. 

We now have that t(z) and r(z) are analytic inside the unit circle, but what 
about on the circle itself and a bit beyond? To show that they are, we need to 
refine our arguments. 

Lemma 7.1. (1) (e~ ie I - P G UP G )\u) = if and only if \u) is an 1? 

eigenstate ofU with support in the edges of G. We call these bound states. 
(2) // |w) = Pq\ — 1,0), then (u\w) — 0, where \u) is any bound state. 

Proof. (1) 

= (e- l6 I-P G UP G )\u) 

(7.9) = e- l0 P^\u)+e- ie P G \u)-P G UP G \u), 

so P G \u) = and the support of \u) lies in the edges of G. Therefore, 
P G U\u) = exp(— i8)\u), which implies that P G U\u) = U\u), and \u) is the 
desired eigenstate. 
(2) This follows from 

(7.10) (u\w) = (u\P G U\ -1,0) = (U*u\ -1,0) = 0. 

□ 

If e~ lB corresponds to a bound state, then the operator (er lB \ — G) is not 
invertible. At the moment this leaves us without a method of defining t(z) and r(z) 
for values of z corresponding to bound states. We would like to show, however, 
that for these values of z, scattering solutions, that is solutions to Eq. I|7.6(l . also 
exist. These solutions can then be used to define t{z) and r(z). 

Let |tti), . . . \uk) be the eigenstates of U, where \uk) e H G and U\uk) = 
exp(— i9k)\uk). Let TL^ be the span of {|iti), . . . \uk )}■ 7~tb is an invariant sub- 
space of both U and G. Let Hi = Tt^~ as a subspace of Tt G , and let Pi be the 
orthogonal projection of H G onto Hi, and Gi = P\GPi. We wish to solve the 
equation 

(7.11) (I Wl - zG 1 )v = zw, 

on Hi- Because of the construction of Hi, we are guaranteed that there exists an 
e > such that (1-^ — zGi) is invertible on Hi for \z\ < 1 + e. This allows us to 
define 

|* G (z)) = v G (z) = {I Hl - zGi^zw 

r(z) = z((0,-l|[/|-l,0) + (0,-l|[/|* G (z))) 

(7.12) t{z) = z((j,j + l\U\-l,0) + (j,j + l\U\V G (z))), 
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and 

-l -l 
\9(z)) = z (k+1) \k,k + l) + r{z)z-^ k+1 ^\k + l,k) 

h— — co A;— — oo 

oo 

(7.13) +|*g(*)> + E *(*)* fc_J >> * + !)■ 

k =J 

We shall summarize our results in a theorem. 

Theorem 7.2. (1) The complex-valued junctions r(z) and t(z), and the 

function \^>q(z)) , which takes values in Hi, are analytic functions on the 
disc \z\ < 1 + e. 

(2) For each z G D = {z\0 < \z\ < 1 + e}, the punctured disc, \^(z)} is an 
element of L°°(E), and \^(z)) is an analytic function of z on the domain 
D with values in L°°(E). 

(3) |^(2)) is the solution to the generalized eigenvalue equation zU\^(z)) = 

!*(*)>■ 

Let us now make two brief observations. First, note that the functions t{0) 
and r(9) are just the restrictions of t(z) and r(z) to the unit circle. Second, we see 
that the above theorem, along with the expressions for t(z) and r{z), provides the 
necessary conditions for all of the integrals in Eq. Ij4.19|l to vanish. 

8. Spectral results 

The formal eigenvectors \"^>(9)), which we just constructed, are elements of 
L°°(E), and hence, for any \v) S L 1 (£'), the inner product (^(6)\v) is well-defined, 
because it is absolutely summable. However, if \v) is merely in L 2 (E), (ty(6)\v) is 
not necessarily summable. Due to the specific form of \^(8)), this sum represents 
a Fourier series with square-summable coefficients, and hence it is a function in 
L 2 (T), where T is the unit circle with the usual measure d9. We shall now exploit 
these two interpretations of our formulas. 

Recall that Tib is the space of bound states, i.e. L 2 eigenfunctions of U. Tib is a 
finite-dimensional, t/-invariant subspace of H.. Each element of Tib has its support 
contained among the edges of G. The function (vector) |^(^)) constructed in the 
last section is seen to be the unique element of L°°(E) such that, 

(1) 17|*(0)> = e-*>|*(0)>, 

(2) (*(0)|-1,O) = 1, 

(3) (*(0)|u) = for any \v) G H b . 

The analytic extension of \^(6)) to D = {z\0 < \z\ < 1 + e} is given by \^(z)} 
(see Eq. (|7.13|n . and a direct application of the Cauchy integral theorem to ^(.z)) 
implies that 

(8.1) ±.j\ 9 \^(e)) = \~i,o). 

If \v) G L 1 , then (^>(9)\v) is a continuous function on the circle, and 

(8.2) ±.J**d0\*(O))(MO)\v) eL°°(E). 
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In particular, if \v) = \ — 1, 0), we have 

(8.3) |-l,0) = i-^ 27r ^|^))(^)|-l,0). 

From this we have that 

U k \-1,0) = -!- / d9e- lk0 \^{e))(^{9)\ - 1,0) 



2tt 
1 

2~/7 



2tt 



d9\^{9))((U^) k ^(9)\- 1,0) 



(i 



271 



(8.4) = hj Q <®W))W)\U k \-l,0). 

Hence, for each \v) in the linear span of the states U k \ — 1,0), where k is an integer, 
and for any trigometric polynomial /, we have the L°° formulas 

(8.5) \v) = i- ( 2 * d6\^{6))^{9)\v) 



2tt 



2?r 



(8.6) f(U)\v) = i- y d0/(e- e )|*(0))(*(e)|i;) 

(8.7) = 4 /^^/(e-^K^Wb)! 2 . 



2tt 



Eq. (|8.6|) is an L°° formula, however, Eq. (|8.7|l is more general. It holds for any 
\v) that is a finite linear combination of the U \ — 1,0), for k an integer. However, 
9 — > is an L 2 function on the circle for each \v) G L 2 {E). Let Hi be 

the closure of the linear span of the U k \ — 1,0) in H — L 2 (E). Hl is a closed 
[/-invariant subspace of H, and we see that Eq. i|8.7|l holds for \v) S Hl and for 
any complex-valued continuous function on T. Thus, by the Riesz-Markov theorem 
there exists a unique measure [i v , a Borel measure on the circle, called the spectral 
measure associated to \v), such that 

(8.8) (v\f(U)\v) = dnv(9)f(e- W ), 

Jo 

and hence, 

(8.9) d , v = friyiy 

27T 

If we define V L :H L -+ L 2 (T) by 

(8.10) Vb|«) = -J=(V(9)\v), 



then Eq. 1|8.7[) implies that ||u|| = ||Vlv||, so that Vl is an isometry from Hl to 
L 2 (T), and Vl{Hl) is a closed subspace of L 2 (T). Finally, we have 

(8.11) {V L UV L - X )g{e) = e-* 9 g(9), 

for g € Vl(Hl) which implies that Vl(Hl) is an invariant subspace of the multipli- 
cation operator g{9) — > exp(—i9)g(9). Therefore, we have that Vl{Hl) — L 2 {T), 
and hence Vj, is a unitary operator of onto L 2 (T). We can summarize our 
results to date in the following theorem. 
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Theorem 8.1. LetHL be the closed, linear, U-invariant subspace ofH = L 2 (E) 
generated by the states U \ — 1,0), for k an integer. Let \^f(6)) be the generalized 
eigenf unctions constructed in the previous section, then for each \v) <G Hl, 

(8.12) V L \v) = J={*(9)\v), 



defines a unitary operator from Hl onto L 2 (T), such that for each g{9) e L 2 (9), 

(8.13) {VlUV^WB) = e- ie g{6). 

Hence, we have constructed part of the spectral decomposition of U. 

Let Hr be the closed, [/-invariant subspace of H generated by the vectors 
U k \j + l,j), for k an integer. In the same way as was done in the previous section, 
we can construct generalized eigenfunctions 

-l 

l*r(0)> = E t r {e)e-^ e \k + i,k) + \^ Gr ) 

k— — oo 
oo 

(8.14) J2( e ~ i(k ~ i)e \ k + !> fc ) + r r (9)e l{k -^ e \k, k + 1)), 

k=j 

which satisfy the results of the last theorem of the previous section. If we follow 
the arguments above, we arrive at the following theorem. 

Theorem 8.2. Let Hr and \^ r (9)) be as above. Then for each \v) G Hr, 

(8.15) V R \v) = -j={* r (6)\v), 



defines a unitary operator from Hr onto L 2 (T), such that for each g(9) e L 2 {9), 

(8.16) (V R UV^)g(e) = e- ie g{6). 

It has already been shown and is easy to see that Hl, 7~Lr, and Hi, are mutually 
orthogonal, [/-invariant subspaces of H = L 2 (E). The final result of this section is 
that these subspaces exhuast H, and, therefore, this completes the spectral decom- 
position of U on H. 

Theorem 8.3. H = L 2 (E) = Hl © Hr Hb ■ 

PROOF. There are many ways to prove this. The following is in the spirit of our 
original argument. If (Hl © Hr © Hb) ± = Ho ^ 0, then we can find a generalized 
eigenfunction |^(^)), such that 

(8.17) U\*(0)) =e- t6 \V(9)), 

where |^(^)) has no overlap with | — 1, 0), + or Hb- Then a direct calculation 
shows that \^(9)} viewed as a function defined as on the edges, E, is identically 
zero. □ 

9. Conclusion 

We have explored a number of properties of quantum walks on graphs. The 
main result is that there is a connection between the S matrix of a graph and the 
properties of a quantum walk on that graph. In particular, the probability that a 
quantum walk traverses a graph in n steps and the conditional hitting time to make 
a walk traversing a graph can be expressed in terms of the transmission amplitude 
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of that graph. This required studying the analyticity properties of the reflection and 
transmission coefficients of the graph. In addition, we have shown that graphs can 
have bound states that quantum walks starting on an edge leading into the graph 
cannot access. We also defined what it means for a quantum walks on a graph 
to be time-reversal invariant, and used this property to show when transmission 
amplitudes for particles coming from different directions will be the same. 

It should be possible to extend all of the arguments here to the case in which a 
graph has more than two tails attached to it. That means that it can have multiple 
entrances and exits. This will simply result in the S-matrix of the graph having 
multiple channels. 



Appendix A 

We want to solve the equation U\^} = e~ t9 \^). Substituting in the explicit 
expression for |<]/), we find, by looking at the indicated edges 

I '- 1 ) 

-if? 2/1,1 \ 

3O0A 



(9.1) 



|0,L4) 
\1A,Q) 
|L4,2) 
|2,L4) 
|2,3) 
|0,1B) 
|lfl,0) 
2) 

|2,li?) 



e- t6 b-i = -i + |(fo A + feos) 
e- ie a 0A = |(1 + 60s) - 3-&0 
e~ te b 0A = hA 



e~ ie a 1A 



aoA 



\aiA 



2 r, 



e~ t6 a 2 = \{a\A + a 1B ) 



e~ l9 a 0B = §(1 + & a) - ^b Q 
e- te boB = biBe* 
e-' l9 aiB = aose^ 
e- ie b 1B 2 - 



3CHA — jo>ib- 

Our object is to find 02 and &_i from which we can find the transmission and 
refelection coefficients, respectively, of the graph. Keeping this in mind, these ten 
equations can readily be reduced to six 



«0B 



e lf> a A 



e b Q A 



h)A) - -^b QB 



(9-2) 



-3»e- 



'OB 



1 2 

-3 + ^{boA + boB) 

2 l 
3(1 + b B) - 3&OA 

1 2 7/> 

-gftOA + ^ e a 0B 

2 • / 

-(a A + e^a 0B ) 



-a 0A - -e v a 0B - 



These equations can then be solved, giving us 



«2 



Ae~ w {\ + e i(t> ){\ 



-i(4e+<t>) 



-Aid 



(l + e"^) 2 - (3e- J ( 40 +« - l) 2 



(9-3) 



'(1 



') 2 +e^(3e- 



-1(49+4,) _ 



l)(e 



3) 



e -U6(i + e -»0)2 _ ( 3e - t (4e+0) _ ^2 
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Note that ai is just the transmission amplitude, t{6), for the graph, and 6-1 is its 
reflection amplitude. 

Appendix B 

We want to show that Eq. (|4.17() is true. A graph is specified by its set of 
vertices, V, and its set of edges, E. A state € L 2 (E) can be expressed as 

(9-4) l^>=E E C (B,A)\B,A). 

Consider a set of vertices in a graph, S, and its complement S. Defining 

p sW = E E hBM 2 

(9-5) P-sW = E E \<\b,A)\ 2 , 

AeS {B,A)£lj a 

we have that 

(9-6) U\\ 2 = PsW+PsW- 

Note that if H^H 2 = 1, then Psiip) is j us t the probability that a particle in the 
state|V>) will be found on one of the edges in the set {(B, A)\A £ S, (B, A) £ uis}, 
with similar considerations for P-g, and this justifies the notation. 
We can also define 

p sW = E E hA,B)\ 2 

Aes (A,B)er A 

(9-7) P-sW = E E \ c (a,b)\ 2 , 

AeS (A,B)Eta 

and we then have that 

(9.8) Uf^PsW+PsW- 
Because the operator, U , is unitary, we clearly have that 

(9.9) P s (j>) + Pg(il>) = P s m) + PgiU^). 

We can now apply this to the state |v&g) -I- | — 1 , 0} and choose the set S to be 
the vertices of G. Application of U to this state gives us 

tf(|* G ) + |-l,0» = e- l9 (|* G )+r(0)|O,-l) 

(9.10) +t(8)\j,j + l)), 

and a short calculation shows that Pg(ip) = Pg(Uip) = 0. We also find that 

PsM = ||* G (0)ll 2 + i 

(9.11) Psm) = ||vI/ G (0)|| 2 + |r(0)| 2 + |^)| 2 . 

Inserting these results into Eq. I|9.9|l Eq. (|4.17|) follows. 

We can actually do a bit more than this. Consider Eq. (|6.7(l in which the 
eigenstate for a particle coming from the left and the eigenstate for a particle 
coming from the right are given. The result in the previous paragraph is just 

(9-12) |n(0)| 2 + M0)| 2 = l, 
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and application of Eq. (|9.9(l to the state \^Gr) + \ j + 1> j) gives us 
(9.13) \r r (8)\ 2 + \t r (9)\ 2 = 1. 

In addition, if we have the state 



(9-14) l^> = E E d {B , A) \B,A), 



we can define 



E E d (B,A) C (B,A) 



AeS (B,A)£u A 



(9.15) 



E E d (A,B) C (A,B)- 



AGS (A,S)er A 



The fact that U is unitary implies that 



(9.16) I s {cf>,i,) + Is{4>,ij) = i s {U<t>,Ui ) ) + I-s{U<t>,U^). 



Applying this equation to the states = \&gi) H- | — 1 , 0} and tp = \&Gr) + \j + 
gives us that 



which relates the incoming amplitudes to the outgoing amplitudes, is a unitary 
matrix. 
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